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Adiabatic approximations are a powerful tool for simplifying nonlinear quantum dynamics, and 
are applicable whenever a system exhibits a hierarchy of time scales. Current interest in small 
nonlinear quantum systems, such as few-mode Bose-Hubbard models, warrants further development 
of adiabatic methods in the particular context of these models. Here we extend our recent work 
on a simple four-mode Bose-Hubbard model with two distinct dynamical time scales, in which 
we showed that among the perturbations around excited stationary states of the system is a slow 
collective excitation that is not present in the Bogoliubov spectrum. We characterized this mode as 
a resonant energy exchange with its frequency shifted by nonlinear effects, and referred to it as a 
second Josephson oscillation, in analogy with the second sound mode of liquid helium II. We now 
generalize our previous theory beyond the mean field regime, and construct a general Bogoliubov 
free quasiparticle theory that explicitly respects the system's adiabatic invariant as well the exact 
conservation of particles. We compare this theory to the numerically exact quantum energy spectrum 
with up to forty particles, and find good agreement over a significant range of parameter space. 

PACS numbers: 03.75.Kk, 03.75.Lm, 67.25.dt 



I. INTRODUCTION 

Interest in the nonlinear quantum dynamics of small 
systems has grown in recent years, as theory and exper- 
iments together approach the mesoscopic regime. With 
Bose-Einstein condensates (BEC) in optical lattices now 
available in many laboratories, few-mode Bose-Hubbard 
(BH) systems provide a realization that is both clean 
and rich [I-JiJ. From a fundamental point of view, small, 
isolated quantum systems have long served as testbeds 
for studying equilibration and thermalization [TUl414j . 
Recent theoretical work has sought to reformulate the 
basic concepts of equilibrium [T51 - [T5] and equilibration 
[TOl [20] • The essential challenge remains great, how- 
ever: Even for quite small quantum systems, the com- 
plete energy spectra can be computationally unattain- 
able. Such standard approximations as Bogoliubov lin- 
earization normally only yield good results in the neigh- 
borhood of the ground state, but a full understanding of 
thermal effects must encompass excited states as well. 

Adiabatic methods may provide an additional theoret- 
ical tool for understanding BH systems, especially since 
it is experimentally possible to tune tunneling rates over 
large ranges in order to impose a time scale hierarchy. In 
our previous work j2lj we introduced a system of two of 
such two-mode BH subsystems coupled as to represent 
a four-mode BH system with a large time-scale separa- 
tion. Each subsystem in this model is a pair of bosonic 
modes coupled (for example, by tunneling between two 
wells) and thus represents an idealized Josephson junc- 
tion [22 -26 . The time-scale separation is implemented 
by a clear hierarchy of the tunneling frequencies coupling 
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the four modes, i.e. the intra subsystem coupling 51 is 
chosen to be by far larger than the inter subsystem cou- 
pling uj; cf Fig. [T^. Therefore the subsystem character- 




a) b) 

FIG. 1. (Color online) Two ways of looking at the same four- 
mode Bose-Hubbard system: a) A'' conserved atoms occupying 
four modes perform multiple Josephson oscillations with a dis- 
tinct time scale hierarchy cj/fl 1; b) In an adiabatic sense, 
the total number J of high-frequency Josephson excitations 
('josons') is also conserved. Hence we may also regard the full 
system as consisting of A'' atoms and J josons, all interacting 
nonlinear ly within two weakly-coupled subsystems. 

istic timescale - the Josephson frequency for each of the 
more tightly coupled mode pairs - is to be the short- 
est one in the entire system. In [3T] we showed that, 
while the total particle number TV is exactly conserved, 
in adiabatic theory there also exists another conserved 
quantity, the total number J of Bogoliubov excitations 
of both of the two-mode BH subsystems. Since these are 
the elementary excitations of idealized Josephson junc- 
tions, we referred to them as 'josons'. Just like atoms, 
josons can be distributed arbitrarily between the two sub- 
systems, but their total number is conserved [see Fig.fTb]. 
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Thus the subsystems may exchange both atoms, by or- 
dinary ('first') Josephson oscillations, and josons (i.e., 
high-frequency energy) , in a second low frequency collec- 
tive mode that we termed 'second Josephson oscillations', 
by analogy with first and second sound in helium II. 

Whereas in [21] we confirmed our adiabatic results with 
comparisons to numerical Gross-Pitaevskii mean field 
theory, in this paper we now explicitly address the quan- 
tum energy spectrum of our four-mode system, and show 
how even highly excited levels can be given accurately by 
our generalization of standard number-conserving Bogoli- 
ubov theory. We then drop our previous restriction to 
the mean field regime of large and relax the restric- 
tion J <^ N. In the spirit of our previous work we are 
then able to derive an extended Bogoliubov theory which 
treats atoms and joson quasiparticles similarly. Imple- 
mented via Holstein-Primakoff transformations (HPT) 
[27] . this formalism ensures the conservation of the to- 
tal number of atoms N as well as the total number of 
josons J, and is thus both technically and qualitatively 
somewhat different from what is usually called number- 
conserving Bogoliubov theory j28 - 31j. Standard Bogoli- 
ubov theory only works well in the neighborhood of the 
^''-particle ground state of the system, but our nonlin- 
early resummed version covers the neighborhood of the 
ground state of each energy band for given N and J. 
This validity range includes many highly excited states 
of the full four-mode system. 

The extensions we will report here provide a useful 
counterpart to other recent work on this same rich four- 
mode system. In |32j . a truncated Wigner approach was 
used for dynamical calculations. Although the so-called 
truncated Wigner approximation only represents non- 
classical dynamics in a limited sense, since it uses strictly 
classical (Liouville) evolution for the Wigner function, 
the results in (32j should be accurate in the mean field 
regime, and otherwise give at least some indication of 
quantum effects. Although [351 attempts to test our pre- 
diction of second Josephson oscillations, the specific cases 
reported there either have J essentially zero, so that sec- 
ond Josephson oscillations do not occur because there 
are no josons present, or else begin with very large dif- 
ferences in the populations of the four modes, so that 
no linear theory can be expected to perform well, and so 
that their Fock state simulations may be sampling many 
trajectories close to unstable classical fixed points, for 
which mean field theory breaks down unusually quickly. 
In this sense [32] and our earlier paper are complemen- 
tary rather than conflicting, and this paper is comple- 
mentary to both, in providing fully quantum mechanical 
results for small perturbations around many excited but 
stationary states. 

This paper will be arranged as follows. In the first 
section we will give a short review of standard Bogoli- 
ubov theory for the simple case of a bosonic two-mode 
system. In the following section we will derive the ex- 
tended version of our theory for larger values of J. The 
third section will be devoted to the structure of the quan- 



tum spectrum and reveals the double ladder structure of 
the energy bands. The last section will focus on the de- 
pendence of the levels on the interactions, present the 
numerical results, and discuss the breakdown of the free 
quasiparticle theory. We will close with a general discus- 
sion. 



II. A REVIEW OF BOGOLIUBOV THEORY IN 
THE TWO-MODE SYSTEM 

The simplest case of the BH model is its two-mode 
version. This simple model has been very extensively 
studied, since it exhibits nontrivial effects such as self- 
trapping, admits a thorough analytic treatment 1221 
f23| [33] and is also experimentally accessible in cold atom 
laboratories [T}]3]. Here it shall serve as a simple and 
well-known example to recall the technique of standard 
(symmetry-breaking) Bogoliubov theory [341 - 137] . The 
two-mode BH Hamiltonian reads as 

H2 = - - (a\a2 + a|ai) + "'^^"'^^ 

1=1 

where Oj , are the annihilation and creation operators of 
the two bosonic modes. A first approximation to such a 
system would be to neglect the interactions completely, 
i.e. to set C/ — >■ 0, and thereby obtain equally spaced 
energy levels (the quanta of Rabi oscillations). However, 
since the interaction energy scales quadratically with par- 
ticle number, this approximation can only be accurate for 
U being very small compared to the coupling fl, which 
however is indeed true in many experimental realizations. 
Nevertheless there exists a much broader regime in which 
the energy levels do remain linearly spaced, to a very 
good approximation; but their spacing can be substan- 
tially different from that given by the simple free particle 
theory with U = 0. This frequency shift depends on the 
interactions and increases with growing U. Bogoliubov 
theory is able to predict these corrections correctly and 
thus rather accurately describes the excitations above the 
ground state even for higher interactions, even though it 
is only a linearization. In Fig.[2]tlie exact energy levels of 
the Hamiltonian ([T]) are compared to the free particle and 
the Bogoliubov approximation. To obtain the Bogoli- 
ubov Hamiltonian one diagonalizes the linear term in ([T]) 
and assumes the symmetric mode to be highly occupied, 
in the sense that its annihilation and creation operators 
can be replaced by the number VN, to yield an effective 
single-mode theory. Keeping only Hamiltonian terms up 
to quadratic order in a and a\ one diagonalizes these 
with a Bogoliubov transformation [571 IMj by introduc- 
ing new bosonic quasiparticle annihilation and creation 
operators b, such that a ~ ub + vW , where — — 1 
to maintain canonical commutation relations. This pro- 
cedure leads to a very simple Hamiltonian of the form 

^Bog = E^ + + UN) b^b (2) 
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state, but above a whole class of states, namely the low- 
est energy states of various energy bands. This extended 
Bogoliuov theory is therefore able to cover a far larger 
part of the spectrum. 



UJ 



UN = 0.1 



UN = 1 



FIG. 2. (Color online) Comparison of the different approxima- 
tions to the exact two-mode spectrum for S7 = 1 and — 30: 
solid red line: exact levels; dotted blue line: free particles; 
dashed black line: Bogoliubov free quasiparticles, all in units 

of no, 



where Eq is an irrelevant c- number constant. Here and 
throughout this paper we assume units such that h = 1. 
The Bogoliubov approximation can therefore be under- 
stood as a theory of free quasiparticles, where the Rabi 
frequency fl of the two-mode system is replaced by the 
Josephson frequency Q = + UN) (which of course 

reduces to fl in the limit J7 — 0). The [/-dependence of 
n is the crucial feature for describing the low-energy ex- 
citations accurately in the presence of interactions. Thus 
Bogoliubov theory is a noninteracting quasiparticle the- 
ory that includes nonlinear interactions among the real 
particles in a mean field regime with one strongly occu- 
pied mode. 

There are several ways to improve the naively classical 
treatment of the highly occupied mode, without signif- 
icantly altering Bogoliubov theory otherwise, by using 
the exact number conservation of the system. In this 
paper as well as in our previous work, we have imple- 
mented number conservation rigorously and exactly by 
using Holstein-Primakoff transformations [21] . which is 
slightly different from the standard formalism introduced 
in [MS]- 

It is important to appreciate that Bogoliubov transfor- 
mations are simply a class of linear canonical transforma- 
tions, in the sense of Hamiltonian mechanics, and that 
there is no intrinsic restriction on the kind of excitations 
to which Bogoliubov theory can be applied. In particular 
it is possible to perform multiple Bogoliubov transforma- 
tions in succession, defining a revised set of higher order 
quasiparticles to absorb certain interactions among the 
quasiparticles that were defined by the first Bogoliubov 
transformation. As we shall see in the following para- 
graphs it is possible in this fashion to construct a gen- 
eralized number-conserving Bogoliubov approach that is 
able to handle excitations not only above the ground 



III. EXTENDED JOSON THEORY 

The four-mode BH system we consider consists of two 
coupled two-mode systems and is described by a Hamil- 
tonian of the form 



n 



(3) 



Hn 



i=l 



Hlr = -— (alj^^aiR + a\^a2R. + H. c. 

with on-site interaction U and index a = L, R for the left 
and right subsystems, respectively. A clear time scale hi- 
erarchy is implemented by choosing Q, the high tunneling 
rate within each subsystem pair, to be far larger than the 
low rate uj between the two subsystems, i.e. restricting 
to the limit tj/fl <C 1. However, as we will see in our 
numerical analysis, this is not actually a very stringent 
condition; the adiabatic results remain remarkably ac- 
curate as long as some degree of time scale separation 
is present. It would not require very much optimism to 
replace the ^ I condition with < I. 



A. Uncoupled subystems 

We will begin by considering uncoupled two-mode sub- 
systems and introduce atom- moving operators , aj, via 
an HPT: 



Na - aiaa CLa = \{°^la + «L)(aia " "2a) • (4) 



These new single-index operators commute with Na — 
X]i=i 2 '^\a'^iot fulfill the usual bosonic commutation 
relations [aQ.,aJj] = 1. Note that since Cq, transfers one 
of the Na atoms, only states with 



a^CLa ^ Na 



(5) 



(i.e. only superpositions of eigenstates of these operators 
whose eigenvalues satisfy this inequality) can be physi- 
cal. With these definitions the Hamiltonians of the sub- 
systems assume the form 

Ha = -^Na + ^alaa + ^Na[Na-2) (6) 



u 



+ ^ ( fllV ~ alcLa + \/ Na~ aaCLaCLa 



At A 



It will be observed that this Hamiltonian provides no cou- 
pling between the physical and unphysical Hilbert spaces. 
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and so the imphcit introduction of unphysical states with 
(aQfla) > (Na) has done no actual harm. 

So far we have followed exactly the procedure of pT] . 
To now go beyond our previously assumed limit J/N ^ 
1, in this paper we must take some care in treating the 
last term in (|6|. In our previous approach we expanded 
this term in inverse powers of \/Na and omitted terms 
of the order of UN~^. But this in fact is a rather bold 
approximation, as will become clear if we do all the trans- 
formations at once. With the help of a second HPT we 
introduce a number-conserving operator a transferring 
atoms between the L and R subsystems, which satisfies 
[a, at] — 1 and commutes with the total particle num- 
ber N — N]^ + Nji and the Oq, modes. For the particle 
number on either side we then get the expression 



L.R 



N ± {a) Vn - ata + - ata a)] . (7) 



On the modes we at first perform a Bogoliubov trans- 
formation [371 [31] letting Ga = uba + vw^ with [ba, b^a\ = 
1, where we leave u and v arbitrary at the moment. These 
parameters will be computed later self-consistently. Since 
we demand the adiabaticity constraint w <C and are 
only interested in dynamical frequencies on the order of 
w, we may apply a rotating wave approximation (RWA) 
discarding all terms not commuting with the number of 
Bogoliubov excitations, i.e. josons, 6|j6q, on each side. 
These are simply the elementary excitations of the two 
idealized Josephson junctions. In this regime the total 
joson number 



J = b\bL 



(8) 



is an adiabatic invariant and may be treated as a new 
constant of motion. This motivates us to finally define a 
joson-moving operator 6, by one last HPT, according to 



j-irbb^hbl 



(9) 



completely analogous to Q. As a direct consequence 
of ([5]) we get the restriction J ^ for physical states, 
since josons are carried by atoms after all, as they simply 
reflect different atom distributions. We may then express 
ajjOa in term of the newly deflned modes d and b, yielding 



dicia 



,,2 r 



J ±{y^\l J -b^b+ \l J -Ubb) 



(10) 

From this it is obvious that for larger J the last term of 
([g]) may more accurately be expanded in inverse powers 
of ^/n — J, where 



(11) 



is the scaled joson number. So far, however, the Bogoli- 
ubov parameters have been involved in the transforma- 
tions, but have not been properly defined yet. This has 



to be done now, self-consistently, by solving a third order 
equation for v? -\- , 



n{n + UN) ■ 



U^N^ 



+ + UN) - 



2\3 



~WJ{u'+v') 



U'J' 



WJ 



U^NJ 



{u^+v^)^ (12) 



{u^+v'), 



while keeping the normalization condition v? — = 1. 
For compactness in the following we will continue to write 
u and V as parameters, but one should keep in mind that 
in fact they turn out to be functions of il, C/, N and J . 
Thus, the leading term of the Hamiltonian is just of the 
standard Bogoliubov type, i.e. 

Hl + HR = nj + 0{Un°) (13) 

with the slightly modified Bogoliubov frequency 



n = Jn{n + u{N - J) 



(14) 



which also has to be determined self-consistently. Note 
that it is not identical with the standard Bogoliubov fre- 
quency of the two- mode system in ([2]), but is shifted 
in the presence of josons. Being interested in the low- 
frequency dynamics only, we ignore the constant parts of 
the Hamiltonian and up to fourth order in the da modes 
finally obtain in terms of atom- and joson-moving oper- 
ators 

U 



Hl + Hr^ -{d^ \/N - a^a+ \/N - d^da)^ (15) 



U 



{u + v)^{u'^ +v^)Ci{P\/ J -bib+^ J -b^bb)^ 



+ ^{u^ + v^)Ci{dWN - dtd + \/n - dtdd) 



:{y^ \J J -m + \J J - mi) 



where the constant Ci = 1 - 0{v'^l{N - J)^). It is to be 
pointed out that the last term in ([l5) is an atom-joson 
coupling. This coupling may be considered as a lingering 
reminder of the fact that, although josons and atoms can 
be considered as independent conserved particles in the 
low-frequency regime, fundamentally josons are motional 
excitations that must be carried by atoms. The atom- 
joson coupling is thus a matter of definition as much as 
of dynamics. 



B. Turning on the coupling 

Let us now turn to the coupling term H^r of the to- 
tal Hamiltonian ([s]), and perform the transformations of 
the preceding section once again. Here once more terms 
occur that need careful treatment. (1 — dl^da/Na)^/"^ , 
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for instance, must be expanded in powers of — J, as 
discussed above. To be consistent we expand again up to 
fourth order in the modes. In RWA we then obtain a 
coupHng Hamiltonian of the form 



(16) 



UJ 



2 

UJ 

2N 



{u^ + v^)C3{aWN - ata + \/ N - d^aa) 



[u^ + v'^){a^^/N - d^d- \/N -d^dd) 



+ -(C3 - C2/iV)(aVA^ - a^a + ^/ N ~ d^ddf 



{u'+v'')C^{y^\lJ- 



b^b- 



J-Wbbf 



where we introduced the constants 
C2 = l-0{v''/N), 
C3 = -^^+Oivy{N - 



(17) 



N-J 



J? 



Thus we have computed the whole four-mode Hamiho- 
nian in terms of atom- and joson-moving operators d and 
b. Note that in contrast to our derivation in ^21j there 
appear terms which vanish in the hmit J/N — >■ 0, for 
example those proportional to y/jjN . These terms are 
necessary for ensuring that the whole transformation per- 
formed on the original operators is canonical; if they are 
omitted in the small J limit, then the operators only ful- 
fill the canonical commutation relations up to corrections 
of the order of 1/7V. This observation does not invalidate 
our derivation in |21j . since there we consistently omit- 
ted terms of that order, and hence could also omit the 
sub-leading terms in J /N . By taking these terms into 
account, now we may also treat states with large values 
of J « iV. 

To now identify the collective excitations we may sim- 




(a^-a)(St-6) 



ply linearize ( 16 ) in d and b which yields 
,.jPb+'^''' 
-N{d^ +af + ^J{P +b) 
'NJ{d^ +d){b^ +b) 
where we used the abbreviations 

u}n=ojC2, ujj ^ uj{u'^ + v'^), ujnj = -^ 

Un = U + uj{C3-C2/N) 
Uj = uj{u^ + v^)C3 - 2U{u + vY{u^ + v'^)Ci 
Unj = U{u + vfCi - uj{u^ + v^)C3. 



(18) 



(19) 



The Hamiltonian ([18j) now just represents a coupled lin- 
ear two-mode system which can immediately be diagonal- 



ized by another Bogoliubov transformation of the form 



U 



V 



(20) 



where the matrices with the Bogoliubov parameters have 
to satisfy the standard constraints UU^^ — VV^^ — 1 and 
UV'^ = VU"'" to ensure canonical bosonic commutation 
relations [d' ,d''^] = 1 and = 1 for the new opera- 

tors. Therefore the Hamiltonian (18) can be written in 
the simple form 

-ffiin = w+a'^a' + Co^y^y = Cj+n + j (21) 



with two independent quantum numbers n, j. The decou- 
pled collective modes, i.e. the eigenfrequencies of (18 1, 



may easily be determined using the Heisenberg equations 
of motion id' = [d',H\i-a\ = ui+d' and ib' = [b',H\in] = 
uj-b' respectively and turn out to be 



NJ 




-,1/2 

+ ujIj{loj + UjJ){ujN + UnN) 
Here we introduced the atom Josephson frequency 



U)N = V^Ar(wjv + UnN) 

= ^/uj{uj + UN) + 0{ujU) 



(22) 



(23) 



which is equal to the lowest mode of the standard Bo- 
goliubow formalism up to corrections of the order of ujU ^ 
and the new Josephson frequency for josons 



= \/ujj{ujj + UjJ) 



(24) 



For the rest of the paper we will refer to this linearized ef- 
fective theory (21) as extended Bogoliubov theory, since 



it is a quasiparticle theory that includes interactions 
among the josons in a similar way as standard Bogoli- 
ubov theory includes interactions among the atoms. 



IV. THE SPECTRUM 



Since the coupling terms in ( 18 1 are rather small, the 



mode- mixing due to the Bogoliubov transformation ( 20 ) 



is also small, and to a good approximation n counts the 
atomic excitations, while j counts the joson excitations. 
Thus, in this adiabatic linear limit the spectrum of the 
full four-mode quantum system is organized as follows. 
The energies are grouped into bands with fixed quantum 
numbers N and J. Since the joson excitations result from 
an unequal distribution of the atoms over the modes, not 
all (n, j) combinations are physical, because J ^ N has 
to be fulfilled. Each of the bands can be described by 



a Hamiltonian of the form (21) and therefore assumes a 
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'pyramidal' structure due to the range of allowed values 
of n and j . 

To get the numerical results for the full four-mode 
quantum systems, we numerically diagonalizc the Hamil- 
tonian. The dimension of the Hilbert space is dim('H) — 
{N + 3) ( + 2) {N + l)/6; desktop computing power has 
allowed us to handle total particle numbers up to 40. 
The pyramidal band structure predicted by the adiabatic 
approximation is immediately recognizable as a close ap- 
proximation to the exact spectrum, to the point where 
we can unambiguously label each exact eigenstate with 
the quantum numbers (iV, J, n, j) of its adiabatic approx- 
imation. The energy eigenstates are of course all simul- 
taneously eigenstates of N. The adiabatic invariance of 
J is reflected in their having expectation values of J that 
are all very close to integers, with very small quantum 
uncertainties. For each exact energy eigenstate, there- 
fore, we can identify its {N, J) band according to the 
adiabatic theory by defining the exact state's J quantum 
number to be the rounding of the exact expectation value 
( J) to the nearest integer. The states within each band 
can be assigned their internal band quantum numbers 
\n,j) in similar fashion, by rounding the nearly-integral 
expectation values of the corresponding operators. 

To quantitatively analyze the mode frequencies, let us 
now turn to the spectrum of the full four-mode system 
itself. This is in principle exactly similar to the analy- 
sis for the two-mode case presented in Sec. We may 
now compare the numerically exact spectrum to the free 
quasiparticle spectrum from the standard Bogoliubov ap- 
proach, and our extended Bogoliubov theory, which in- 
cludes interactions among the Bogoliubov quasiparticles. 
Since for fixed N the spectrum is organized in J-bands, 
we may restrict ourselves to one of these, and group 
the states according to their n, j-excitations. The re- 
sult for N = 10, ( J) = 3 and the moderate interaction of 
UN = 0.05 is shown in Fig. [3] The 'pyramidal' structure 
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FIG. 3. (Color online) Comparison of the different approxi- 
mations to the exact four-mode spectrum (for A'^ = 10, J = 3, 
Q = 1, (jj = 0.1 and UN — 0.05): solid red line: exact lev- 
els; dashed black line: extended Bogoliubov; dotted blue line: 
standard Bogoliubov free quasiparticles, all in units of HQ 



which appears correctly resembles the linear Hamiltonian 
with two quantum numbers. This structure is truncated, 
however, since for N, J fixed both n and j are bounded 
from above, as they must fulfill ^ n ^ N — J and 
^ j ^ J, respectively. (This restriction is of no prac- 
tical consequence, since the extended Bogoliubov theory 
breaks down for these highly excited states anyway. Nev- 
ertheless it is possible to derive a complementary theory, 
quite analogously to this one, which correctly describes 
the upper end of each J-band, leaving poorly described 
only the middle ranges of large bands.) It can be seen 
that for low excitations within any band the extended 
Bogoliubov theory is a clear improvement over the stan- 
dard Bogoliubov approximation. 



V. FOCUS ON INTERACTIONS 

Let us now turn to the dependence of the energy lev- 
els on the interaction strength. To compare quantum 
rather than mean field effects for different total num- 
bers of particles it is convenient to keep the product U N 
constant. We may now numerically calculate the energy 
levels of each J band and compare them with standard 
Bogoliubov and our extended Bogoliubov theory. Typical 
results are shown in Fig. |4] Part a) shows the numeri- 
cally exact levels as solid red curves and the extended 
Bogoliubov theory in dashed black curves. The lowest 
energy level of the band has been set to zero for clar- 
ity. For moderate interaction strengths the agreement 
between the exact levels and the extended Bogoliubov 
theory is obvious, even up to higher excited states in the 
band. For even higher excitations, linearization clearly 
fails, and interactions between the quasiparticles have to 
be taken into account. In part b) the exact levels (again 
solid red) are compared to standard Bogoliubov (dashed 
black) . The exact lower frequency spacing in the pyramid 
is evidently shifted below the Bogoliubov beat, as we have 
argued, and only approaches the Bogoliubov value in the 
limit C/ — > 0. For the higher of the two pyramid spacings, 
conversely, the standard Bogoliubov theory consistently 
gives a substantial underestimate. The extended Bogoli- 
ubov theory is in contrast able to capture the main fea- 
tures of the spectrum very well, despite being a simple 
linear approximation with only a two-frequency pyramid 
spectrum. In that sense our model is linear but non- 
trivial, since it includes more effects of the interactions 
in a resummed way via the nonlinear transformations. 

While the numerical problem becomes very hard for 
larger particle numbers than those presented here, ex- 
perience with Bose-Hubbard systems generally indicates 
quite good quantum-classical correspondence, at least 
within the superfluid phase, so that the numerical Gross- 
Pitaevskii simulations presented in [51] should accurately 
represent the quantum dynamics for larger N. The good 
behavior of our theory in that regime, together with the 
successes shown here, therefore confirm that it is valid 
over a wide range of parameter space. Like any simple 
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FIG. 4. (Color online) Exact spectrum of the four-mode sys- 
tem for A'^ = 36, fl = 1 and uj = 0.1, shown is the J = 18 
band; a) solid red curve: exact diagonalization, dashed black 
curve: extended Bogoliubov theory; b) solid red curve: exact 
diagonalization, dashed black curve: standard Bogoliubov, all 
in units of Hfl 



theory, however, it does have its limitations. 

As we pointed out in our previous paper, the inter- 
actions of the josons tend to have the opposite sign to 
the atomic interactions. When atoms repel each other, 
josons attract. This leads to the possibility of dynamical 
instability, and therefore to a breakdown of our theory, 
when the joson frequency turns imaginary at a critical 
point. This opposite interaction of josons and atoms 
holds strictly in the J N limit, but also survives in 
the more general case presented in this paper, since for 
relatively small repulsive interactions 



UN > 



2{u + v) 



0{J/N) 



(25) 



the joson interaction Uj becomes negative according to 
( 19 ). Thus, also in this framework the joson frequency fi- 
nally turns imaginary for high interactions and the model 
breaks down. In particular this means that the number- 
locked regime will not be accessible with this theory. We 
can see moreover that the gradual breakdown of the lin- 
ear joson theory is accompanied by a series of crossings 



among the levels of each band. In particular even the 
lowest two states ultimately approach each other, and 
thus the low frequency goes to zero. At this point the 
level structure no longer consists of two linear modes, 
and is not be expected to be covered by a simple linear 
theory of this kind. 

Another possible limitation to be explored is the re- 
liance on the clear time scale separation. The whole sys- 
tem was constructed in the first place to implement a 
clear hirarcy of the tunneling frequencies according to 
w/n <C 1, to enable the adiabatic RWA to hold. We 
may now also numerically check the validity of the the- 
ory for larger values of the frequency ratio. The limit 
w — > is rather uninteresting, since in the regime of 
moderate interactions the linear part of the Hamiltonian 
then dominates and the whole spectrum is almost linear. 
The limit — > w on the contrary is more interesting. In 
Fig.[5]the J — 1& band is shown for = 32 particles and 

= 0.4 and u = 0.1. In this case the theory still yields 
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FIG. 5. (Color online) Exact spectrum of the four-mode sys- 
tem for N = 32, Q, = 0.4 and to = 0.1, shown is the J = 16 
band; a) solid red curve: exact diagonalization, dashed black 
curve: extended Bogoliubov theory; b) solid red curve: exact 
diagonalization, dashed black curve: standard Bogoliubov, all 
in units of 2.5hi} 



good agreement with the exact spectrum. For smaller 
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values offlJ starts to be not a properly defined quantum 
number since fluctuations rise, which make it numerically 
challenging to extract the band structure correctly, espe- 
cially since the bands overlap significantly in that regime. 
However, we can argue that the adiabatic theory does not 
rely on the strong adiabatic limit, but rather still holds 
for larger values of the frequency ratio oj/Q. So in fact 
the regime of first and second Josehpson oscillations in 
this four-mode system is not a narrow region in parame- 
ter space, but actually the generic behavior in the regime 
of moderate interactions. 



VI. DISCUSSION 

In this paper we have demonstrated a powerful exten- 
sion of the number-conserving Bogoliubov linearization 
theory which can be applied to bosonic many-body sys- 
tems with time scale hierarchies. This resummed lin- 
earization conserves the number of high-frequency exci- 
tation quanta, which in exact terms is adiabatic invariant 
rather than strictly conserved, as well as the exactly con- 
served particle number. While we have demonstrated 
it for a very simple quantum many-body system, the 
general idea should be quite robust, as adiabatic meth- 
ods usually are, and be extensible to more complicated 
systems of many kinds. The main value of the method 
is in providing simple quasiparticle descriptions valid in 
the Hilbert space neighborhoods of many excited states, 
rather than only of the ground state. 

As well as thus demonstrating a calculational tech- 
nique, we consider our study of this four-mode Bose- 
Hubbard system as providing a conceptually instructive 
toy model for an important general phenomenon: the 
quantum dynamics of hydrodynamic collective excita- 
tions. 

Hydrodynamic collective excitations of many-body 
systems close to equilibrium but far from the ground 
state, such as sound waves in air, are usually described 
in terms of classical fields, such as those of temperatine, 
velocity, and pressure. These thermodynamic quanti- 
ties are in fact directly related to microscopic dynamical 
quantities, such as the densities of particles, momentum, 
or energy. The reason why these particular quantities 
play the main roles in low-energy effective theory is that 
they are conserved, and so they cannot be locally changed 
even by very rapid microphysics, but only by the slower 
processes of transport. The description of gas dynamics 
in terms of classical hydrodynamic collective variables 
such as pressure and temperature is therefore ultimately 
justified by a dynamical time scale hierarchy. 

But time scale hierarchies can be exploited very effec- 
tively within pure-state Hamiltonian quantum mechanics 
by using adiabatic theory, without assuming either sta- 
tistical equilibrium or the classical limit. 

This raises the question of whether a purely quantum 
and Hamiltonian theory could be applied in the meso- 
scopic regime. 



Ordinary sound waves in air are, after all, rather sim- 
ple, linear excitations. Might they not have a simple 
quantum description, with canonical kets expressed in 
a basis of quantized quasiparticle excitations, and not 
just as a phenomenological requantization of the classi- 
cal field theory, but rigorously derived from many-body 
first principles? This question may be entirely academic 
for macroscopic systems with a huge number of degrees 
of freedom, but is important in the mesoscopic regime, 
where such theories might be derived directly from the 
many-body Hamiltonian, as in the simple example of our 
four-mode BH system. Thus, this issue is very impor- 
tant for understanding the mesoscopic interface between 
quantum mechanics and thermodynamics, since here the 
new approach can directly be compared to the standard 
thermodynamical treatment. And of course the inter- 
face between quantum mechanics and thermodynamics in 
general is an important frontier of fundamental physics, 
which with the advent of nanotechnology and quantum 
information processing may even turn out to be of prac- 
tical significance. 

The quantum theory of hydrodynamics can be ap- 
proached by considering a volume of gas to consist of 
many adjacent sub- volumes, such that the free motion of 
atoms across the arbitrarily defined sub-volume bound- 
aries provides a coupling between the sub-vohimes. Mo- 
tion and collisions within each sub- volume are then con- 
sidered as internal dynamics for the sub- volume, consid- 
ered as a dynamical sub-system of the whole gas. By def- 
inition, then, each sub-system's internal dynamics con- 
serves the sub-system's total energy, momentum, and 
particle number. If the subvolumes are large enough, 
the rate of proportional change in these quantities, due 
to particles moving between sub- volumes, must be low. 
The effective low-frequency description of the problem 
therefore takes these slowly changing quantities as its el- 
ementary degrees of freedom. 

This means, in particular, that the local sub-volume 
energy is effectively treated in much the same way as 
the sub-volume particle number. The local energy might 
almost be a second kind of fluid, mixed together with the 
fluid of particles. Most of this energy is usually thermal, 
and in this sense the kinetic theory of heat reduces, as 
a hydrodynamic theory at low frequencies, to something 
very much like the previous theory which it replaced - 
the caloric theory of heat developed in the 18th century 
by Lavoisier and Laplace, according to which heat was 
considered a material fluid analogous to electric charge. 

The effective behavior of energy as a dynamical degree 
of freedom is in fact familiar even in dynamical systems 
much simpler than an interacting gas. To observe a pair 
of weakly coupled pendulums in underdamped motion, 
and see how the swing amplitude oscillates slowly back 
and forth from one pendulum to the other, is to notice 
that the energy distribution between the two pendulums 
behaves very much like a pendulum itself. 

The contribution we have made in this paper is to 
study a quantum many-body system that is almost as 
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simple as those pendulums, and show how the low- 
frequency description of local energy as an effective sec- 
ond fluid can be implemented quantum mechanically, en- 
tirely within a canonical Hamiltonian closed system. In 
effect we have simplified the hydrodynamic construction 
of sub- volumes to the minimal case of just two neighbor- 
ing sub- volumes, and enforced relatively slow motion be- 
tween them by requiring it to proceed by slow tunneling. 
We have then further simplified each sub- volume by ide- 
alizing them as two-site bosonic lattice gases, modeling 
their internal sub- volume dynamics in quantum Joseph- 
son form. As a toy model for the hydrodynamic problem, 
then, our system is clearly crude to the point of absolute 
minimality, but does manage to capture some qualita- 
tively crucial features of the real phenomenon. Because 
we have gone beyond the trivial case of non-interacting 
particles, moreover, we have identified some non-trivial 
interaction effects. 

In our previous paper on this system, we showed that 
the {quasi-)particles of the 'heat fluid,' which we have 
called 'josons', tend to attract each other if the funda- 
mental particles repel one another, and vice versa. We 
then showed in consequence of this that for repulsively 
interacting atoms, the attractive interaction of josons 
makes the transport of 'heat' between our two subsystems 
slower than one would anticipate from a naive application 
of the coupled pendulum picture. We have now shown 
how this picture persists into the fully quantum regime. 
In particular we have demonstrated that an extended Bo- 
goliubov theory, based on treating the joson number J on 
the same footing as the atom number N , can yield good 
approximations for the exact many-body energy levels 



within each {N^ J) band. There is a somewhat dizzying 
but pleasing symmetry in the way that the successive Bo- 
goliubov transformations absorb subtler forms of inter- 
action into new forms on non-interacting quasi-particle. 
Just as the first Bogoliubov transformation absorbs non- 
linear corrections to Rabi oscillations into the frequency- 
renormalized Josephson oscillations, so the second Bo- 
goliubov transformation absorbs nonlinear corrections to 
the linear Bogoliubov theory into frequency-renormalized 
second Josephson oscillations. This is, after all, what hy- 
drodynamics is all about: the emergence of simple col- 
lective modes from highly nonlinear microphysics. 

While fully hydrodynamic and thermodynamic phe- 
nomena such as irreversibility and turbulence are obvi- 
ously far beyond the scope of our minimal toy model, 
nonetheless we can already see what may perhaps be 
some first, faint indications of how they emerge from 
Hamiltonian closed-system quantum mechanics. Our nu- 
merical plots of exact energy levels show a number of 
avoided crossings, which our linear extended Bogoliubov 
theory naturally does not recognize, even where it other- 
wise follows the exact level curves very well. In systems 
with many more degrees of freedom, one can anticipate 
that such avoided crossings may perhaps proliferate until 
the free-quasiparticle theory is thoroughly eroded. Bose- 
Hubbard systems such as ours tend in general to be dy- 
namically chaotic, and the application of adiabatic theory 
such as ours in the presence of chaos presents numerous 
difficulties [381I39]. With further research, the nonlinear 
quantum dynamics of simple many-body systems with 
time scale hierarchies may offer further insights into im- 
portant basic problems. 
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